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Abstract: We present the implementation of the dipole subtraction formalism for the 
real radiation contributions to any next-to-leading order QCD process in the Mad- 
Graph/MadEvent framework. Both massless and massive dipoles are considered. Starting 
from a specific (n + l)-particle process the package provides a Fortran code for all possible 
dipoles to all Born processes that constitute the subtraction term to the (n+l)-particle pro- 
cess. The output files are given in the usual "MadGraph StandAlone" style using helicity 
amplitudes. 
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1. Introduction 

Physics studies at the upcoming CERN LHC collider will frequently involve multi-particle 
final states. Especially searches for physics beyond the standard model rely on the recon- 
struction of new particles from their decay products, often through decay chains. Equally, 
requiring accompanying particles in the final state may serve to improve the ratio of signal 
to background processes, as done for example in the Higgs boson search through the vector 
boson fusion channel. Meaningful searches for these signals require not only a very good 
anticipation of the expected signal, but also of all standard model backgrounds yielding 
identical final state signatures. From the theoretical point of view, high precision implies 
that one has to go beyond the leading order in perturbation theory to be able to keep up 
with the precision of the measurements. 

For leading order processes there have been many developments concerning event gener- 
ation and simulation tools in the last two decades such as MadGraph/MadEvent [1-3] 
CompHEP/CalcHEP [4]/ [5], SHERPA [6] and WHIZARD [7] and also programs using 
different approaches such as ALPGEN [8] and HELAC [9]. All these programs are multi- 
purpose event generator tools, which are able to compute any process (up to technical 
restrictions in the multiplicity) within the standard model, or within alternative theories 
specified by their interaction Lagrangian or Feynman rules. They usually provide event 
information which can be interfaced into parton shower, hadronization and/or detector 



Next-to-leading order (NLO) calculations are at present performed on a process-by-process 
basis. The widely-used programs MCFM [10,11], NLOJET++ [12], MC@NLO [13,14] and 
programs based on the POWHEG method [15-20] collect a variety of different processes 
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in a standardized framework, the latter two methods also match the NLO calculation onto 
a parton shower. 

The NLO QCD corrections to a given process with a n-parton final state receive two types 
of contributions: the one-loop virtual correction to the (2 — * n)-parton scattering process, 
and the real emission correction from all possible (2 —* n + l)-parton scattering processes. 
For the numerical evaluation, one has to be able to compute both types of contributions 
separately. 

The computation of one-loop corrections to multi-particle scattering amplitudes was per- 
formed on a case-by-case basis up-to-now, the calculational complexity increased consider- 
ably with increasing number of external partons. Since only a limited number of one-loop 
integrals can appear [21,22] in the final result, the calculation of one- loop corrections can 
be reformulated as determination of the coefficients of these basis integrals, plus potential 
rational terms. Enormous progress [23-42] has been made in the recent past in the sys- 
tematic determination of the one-loop integral coefficients and rational terms, and steps 
towards fully automated programs for the calculation of one-loop multi-parton amplitudes 
were made with the packages Cut Tools [43], BlackHat [44], Rocket [45] and GOLEM [46]. 
The real emission corrections contain soft and collinear singularities, which become ex- 
plicit only after integration over the appropriate real radiation phase space yielding a hard 
n-parton final state. They are canceled by the singularities from the virtual one-loop contri- 
butions, thus yielding a finite NLO correction. To systematically extract the real radiation 
singularities from arbitrary processes, a variety of methods, based either on phase-space 
slicing [47] or on the introduction of process-independent subtraction terms [48] have been 
proposed. Several different algorithms to derive subtraction terms are available: residue 
subtraction [49], dipole subtraction [50,51] and antenna subtraction [52-55]. 
Especially the dipole subtraction formalism, which provides local subtraction terms for all 
possible initial and final state configurations [50] and allows to account for radiation off 
massive partons [51], is used very widely in NLO calculations. It has recently also been au- 
tomated in the SHERPA framework [56] and the TeVJet framework [57], and most recently 
in the form of independent libraries [58] interfaced to MadGraph. The dipole subtraction 
within the SHERPA framework is not available as a stand-alone tool, while within the 
TeVJet framework, the user needs to provide all the necessary process dependent infor- 
mation. Moreover both approaches have only included massless particles for the dipoles. 
There is still no general tool available which is able to produce the dipole terms for an 
arbitrary process and which can also deal with massive partons. 

In this paper we present MadDipole, an automatic generation of the dipole subtraction 
terms using the MadGraph/MadEvent framework. The results are Fortran subroutines 
which return the squared amplitude for all possible dipole configurations in the usual 
MadGraph style. We describe the construction of the dipole subtraction terms and their 
implementation in Section |2| Results from various checks of the implementation are pro- 
vided Section ||, and instructions on the practical usage of the package are contained in 
Section f|. 
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2. Construction of dipole terms 

The fundamental building blocks of the subtraction terms in the dipole formalism [50,51] 
are dipole splitting functions V^fe, which involve only three partons: emitter i, unresolved 
parton j, spectator k. A dipole splitting function accounts for the collinear limit of j with 
i, and for part of the soft limit of j in between i and k. The dipole factors, which constitute 
the subtraction terms, are obtained by multiplication with reduced matrix elements, where 
partons i, j and k are replaced by recombined pseudo-partons ij, k. The full soft behavior 
is recovered after summing all dipole factors. 

Throughout the whole paper we are using the notation introduced in Refs. [50] and [51]. 
Independent on whether we have initial or final state particles we can write an arbitrary 
dipole in the form 

T>ij,k ~ m {l,...ij,...,k,...,m + l\ 2 lJ V ij:k \l,...ij,...,k,...,m + l) m . (2.1) 

ij 

The amplitude factors (. . . | ('bra') and | . . .) ('ket') on the right hand side are tensors in 
color space. The helicities of the external particles in them are a priori fixed (but can be 
summed over for unpolarized processes), while the helicities of the pseudo-partons have to 
be summed over after contraction with the dipole splitting function. 

These Born-level amplitude factors are provided by the usual MadGraph code. The two 
elements that combine the ket with the bra are the additional color structure ^ %3 and 

ij 

the dipole splitting function Vy ^. 

2.1 Color and helicity management 

For the calculation of the color factors there already exist routines in the MadGraph pro- 
gram. Our intension was to use exactly these routines because this code is very well- 
confirmed and efficient. We have included the additional color operator, • Tj,-, by 
rewriting the internal MadGraph color labelling for the ket-side only. After insertion of 
this color operator the color structure is no longer multiplied by its own complex conju- 
gate and therefore the routine that squares the color needed to be altered, to multiply 
the modified ket by its original complex conjugate. We emphasize that due to the facto- 
rial growth of the color factors MadGraph can not handle more than seven colored particles. 

For the insertion of the splitting function Vj,- ^ several changes with respect to the original 
code are required. One has to keep in mind that in general the splitting function is a tensor 
in helicity space, i.e., 

v ii>fc = {n\Vij,kW) = v{£ fc . 

As MadGraph deals with helicity amplitudes, we have to write the dipole in a slightly 
different way to be able to include the calculation of the splitting function in the code. 
Neglecting the color for a moment we start from the definition of the dipole in fl2.1|) and 
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by inserting a full set of helicity states —g^ = J2x e uW e vW we S e t 

(2.2) 
,k, ...,m + l) m 

= rn(- ■ -\x b V(X b , \ a ) X a \ ■ ■ .}m 

with V(X b , A a ) = e M (A b ) V^ fc e*(A a ) and e^(A) M | . . .) m = x\ ■ ■ )m- 

The polarization vectors are calculated using the HELAS routines [59] already available in 
the MadGraph code. The parts that are diagonal in helicity space are trivial to calculate in 
that sense that one only has to multiply the MadGraph output for the squared amplitude 
for a given helicity combination with the splitting function. To calculate the off-diagonal 
helicity terms, the amplitude for each helicity combination is stored and then combined 
with the according amplitude with opposite helicity. 

For the calculation of the splitting functions and for the remaping of the momenta we use 
modified versions of the routines used in MCFM [10, 11]. 

2.2 Massive particles 

If some of the masses of the external particles are non-zero, in particular for processes 
involving top and/or bottom quarks, there are dipoles for which the unresolved parton is 
massive. In this case the collinear singularities are regulated by the mass of the unresolved 
parton and the unsubtracted matrix element does therefore no longer diverge in these 
collinear limits, but only develops potentially large logarithms. Our code still calculates 
all possible dipoles, also in which the unresolved parton is massive, but puts them in a 
separate subroutine, dipolsumf inite ( . . .), that is not evaluated by default. In the limit 
of large center of momentum energy or, similarly, small external masses, the user can easily 
include the non-divergent dipoles to subtract the associated large logarithms, which can 
then be included analytically through the integrated subtraction terms. In the limit of zero 
external masses we have checked that the results obtained after summing all dipoles are 
the same as obtained by generating the code with massless particles from the start. 

2.3 Phase space restriction 

The calculation of the subtraction terms is only necessary in the vicinity of a soft and/or 
collinear limit. Away from these limits the amplitude is finite and there is in principle no 
need to calculate the computationally heavy subtraction terms. The distinction between 
regions near to a singularity from regions without need for a subtraction can be parame- 
terized by a parameter usually labelled a with a 6 [0, 1], which was introduced in Ref. [60] 
for processes involving partons only in the final state. The case with incoming hadrons, 
i.e., with partons in the initial state, is described in Ref. [12]. 
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T>ij,k ~ m(l, -ij, -, k, m + l| M V^ fe „|1, ...ij, k, m + l) m 

= m (l,...ij,...,k,...,m + (-5^) V^ fc (-#/) slh-ij,- 
= £ 4..|,,^(A^)V^4(Ay(A a ),|...) m 



Using the notation of Ref. [12], the contribution from the subtraction term to the differen- 
tial cross section can be written as 



d°ab= dT ( - n+1 \p a ,p b , Pl ,...,p n + l)—^- 

^{n+l} 



{n+1} 



I ^2 V ijAPa,Pb,Pl, ■■■,Pn+l)F ( j a \p a ,Pb,Pl, ~,Pij,Pk, ~,Pn+l)Q(yij,k < Oi) 

(, pairs ky^i,j 

+ E 



airs k^iij 



1^ij(Pa,Pb,Pl, -,Pn+l)F) > (p a ,Pb,Pl, ~,Pij, ■;Pn+l)@{l ~ X ij>a < a) 



pairs 

i,3 



+(a <-> b) 

+ Y1 [ V k(Pa,Pb,Pl, ...,p„+l)Fj n; \p a ,Pb,Pl, -,Pk, ..,Pn+l)&(Ui < <*) + (d 6) 

i+k 

Pb,Pi, ■■■,Pn+i)Fj n \p a ,Pb,Pi, ...,p n +i)@(vi < a) + (a <->■ 6)1 I . 

(2.3) 

The functions Pj^fc, P^-, P^ 1 and D m ' b are the dipole terms for the various combinations 
for emitter and spectator. J2{n+i} denotes the summation over all possible configurations 
for this (n + l)-particle phase space which is labelled as <ir( n+1 ) and the factor S'{ n+1 } is the 
symmetry factor for identical particles. We have introduced four different a-parameters, 
one for each type of dipoles. In our code they are called alpha_f f , alpha_f i, alpha_if and 
alpha_ii for the final-final, finial-initial, initial-final and initial-initial dipoles, respectively. 
The actual values for these parameters are by default set to unity, corresponding to the 
original formulation of the dipole subtraction method [50,51], but can be changed by the 
user in the file dipolsum.f. 

It has to be kept in mind that the integrated dipole factors, which are to be added with 
the virtual n-parton contribution, will also depend on a. For case of massless partons, the 
a-dependence of the integrated terms is stated in [12, 60] while for massive partons results 
for most cases can be found in [61,62]. 



3. Checks 



The MadDipole package provides a code, check_dip . f , which allows the user to test the 
limits of the (n + l)-particle matrix element and the dipole subtraction terms. This code 
builds up a trajectory of randomly selected phase space point approaching a given soft or 
collinear limit of the (n + l)-parton matrix element and yields the values of matrix element, 
sum of all the dipoles, and their ratio along this trajectory. The result is printed to the 
screen in a small table for which each successive row is closer to the singularity. The ratio 
between matrix element and the sum of the dipoles should go to unity. We have tested 
our code in all possible limits, both for massless as well as massive dipoles and found no 
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Figure 1: Matrix element squared |Mr| 2 (upper plots, solid line) and the subtraction terms D 
(upper plots, dashed/dotted/dot-dashed lines) for the process e + (pi)e~ (P2) — ► Z — ► q{p3)q(p4)g{p5) 
as a function of 1 — xq — {ps-Ps)/ (pi-Pz) and x g = 1 — {pz-Pa) / (pi-Pz) in figures (a) and (b), 
respectively. Also plotted are the ratio D/\Mji\ 2 , the difference |Mr| 2 — D (averaged over 100 
random points per bin) and the maximal difference max(||A/fl| 2 — D\) per bin. The dashed lines 
include the dipoles for each point in phase space, a// = 1, while for the dotted a// = O.f and 
dot-dashed a// = O.Of the phase space for the dipoles has been restricted to the collinear/soft 
regions. 



inconsistencies. Choosing small values for a-parameters, e.g., a = 0.1, improves the com- 
putation time and the convergence of the subtraction procedure. 

To show that the subtraction terms are implemented correctly we provide a couple of ex- 
amples in the form of plots and argue that the cancellation between the matrix element 
squared and the subtraction term is as expected. In the figures [l] and [2| we show the matrix 
element squared and the subtraction term as a function of a variable that represents a soft 
and/or collinear limit of the process specified. For these figures we have binned the x-axis 
(equally sized bins for the logarithmic scale) and generated random points in phase space 
to fill each of the bins with exactly 100 events. In the upper plot, |Mr| 2 and D are the per 
bin averages of the matrix element squared and the subtraction term, respectively. The 
second to upper plot shows the per bin average of the ratio of the matrix element squared 
and the subtraction term, while the third plot from the top shows the per bin average of 
the difference. The lowest plot shows the absolute value of the maximal difference among 
the 100 points in a bin. To show the effects of the phase space restriction for the dipoles, 
see section all the plots are given for a = 1 (dashed lines), a = 0.1 (dotted lines) and 
a = 0.01 (dot-dashed lines). 

In figure [l|(a) we show the matrix element squared and the subtraction term as a function of 
1—Xq, where Xq is the fraction of the energy carried by the anti-quark, xq = S34 ^ 45 = 1 — §f|> 
with Sij = Pi.pj. For this process, e + {p\)e~{p2) — ► Z — ► q{pz)q{pA)g{p^)i there are only 
final-final state dipoles contributing to the subtraction term. The center of mass energy 
is set equal to the Z boson mass = mz- To restrict the discussion to the collinear 
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divergence only, points close to the soft divergence (xq = x q = 1) have been removed by 
forcing x q + x q < 1.5 in the generation of the phase space points. 

From the upper plot it is clear that both the matrix element squared and the subtraction 
term diverge in the collinear limit Xq — ► 1, as 1/xq. The ratio D/\Mr\ 2 goes to 1 and the 
average values of the differences fluctuate close to as can be seen in the second and third 
plots from the top. The numerical fluctuations for small 1 — xq can be completely explained 
by statistical fluctuations. They are of the order of 1% of the maximal difference given 
in the lower plot. As can be expected, the cancellations are not exact, which is shown 
by the lower plot. The maximal difference between |M^| 2 and D rises like l/y'l — Xq, 
which does not lead to a divergent phase space integral, because the integration measure is 
proportional to xg. The small peaks/fluctuations in the region for small Xq are due to the 
fact that we are approaching the other collinear limit, i.e., for which the gluon is collinear 
to the anti-quark x g — ► 1, where the matrix elements squared and the subtraction term 
also diverge. 

In figure ||(b) the same matrix elements and subtraction terms are presented, but as a 
function of the fraction of the energy carried away by the gluon X g — i— X q Xq> The 
limit for which x g goes to zero represents the soft divergence of this process, while the 
collinear divergences for this process are removed by excluding phase space points for 
which xq > (1 + x q )/2 or x q > (1 + Xq)/2. The same conclusion as for figure (j](a) can 
be drawn here: the matrix element squared and the subtraction term diverge in the soft 
limit, their ratio goes to one and the average difference to zero, while the absolute value of 
the maximal difference still rises when approaching the soft limit, but does not lead to a 
divergent phase space integral. 

An example for a collinear limit between final state and initial state particles is given 
in figure ||(a). In this plot the matrix element squared and the subtraction term for the 
process e~ (p{)q{p2) — ► e~ (p3)q(p4)g(ps) are given as a function of the invariant mass of the 
initial state quark and the final state gluon S2s/si2- As this invariant mass goes to zero, 
the matrix element squared and the subtraction term diverge like I/S25, an d their ratio 
goes to one. To remove the other possible divergences a cut on the momentum transfered 
S13/S12 > 0.5 and on the invariant mass of the final state quark and gluon S45/S12 > 0.2 
have been imposed. Like for the final-final state dipoles the average difference goes to a 
constant, as can be seen from the second plot from the bottom, but the dipoles have a 
sizeable constant contribution. Therefore the normalization of the average value for the 
difference |M#| 2 — D depends on the number of dipoles included for the phase space point. 
If all the dipoles are included for all points the difference goes to a smaller constant than if 
we restrict the phase space of the subtraction term to be close to the singularities by setting 
a < 1. Due to this restriction only the dipoles to cancel that divergence are included in the 
subtraction term and therefore give a smaller constant contribution, hence the difference 
|M/j| 2 — D is larger. Also here the maximal left-over difference, the lowest plot, increases 
for small invariant masses but does not lead to a divergent phase space integral. 
In figure |2](b) an example with massive final state particles is shown. The process is tt 
production at a linear collider, e + (pi)e~ (P2) — ► t{p3)i(Pi)9(p5) at 1 TeV center of mass 
energy. The plot shows a behavior very similar to the massless case, fig. EL and the conclu- 
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Figure 2: Matrix element squared |Mr| 2 (upper plots, solid line) and the subtraction term D 
(upper plots, dashed/dotted/dot-dashed lines) for (a) the process e~ {px)q(jpz) — * e~ {pz)q{pi)g{p5) 
as a function of S25/S12 = {p2-Ps) / '(pi-Pa) and (b) the process e + {pi)e~(p 2 ) — > Z — > t(p 3 )i(p4)g(p 5 ) 
as a function of x g = 1 — (P3-P4)/{pi-P2)- Also plotted are the ratio D/\Mr\ 2 , the difference 
\Mr | 2 — D (averaged over 100 random points per bin) and the maximal difference max(| |M^| 2 — D\) 
per bin. The dashed lines include the dipoles for each point in phase space, a = 1, while for the 
dotted a = 0.1 and dot-dashed a = 0.01 the phase space for the dipoles has been restricted to the 
collinear/soft regions. 



sions drawn there apply to this plot as well. 

As a further check we have tested the code extensively against MCFM [10,11]. We have 
generated random points in phase space and compared the subtraction terms calculated 
by MCFM with the subtraction terms calculated by our code. See table |l] for a list of 
processes that have been checked. We observed differences only in the case where dipoles 
were introduced entirely to cancel collinear limits, which can be made independently of the 
spectator particle. In our code all possible dipoles are calculated, which implies a sum over 
all spectator particles. However, if there is only a collinear divergence, i.e., the unresolved 
parton cannot go soft, this sum is redundant and one dipole with the appropriate coefficient 
is enough to cancel the singularity. In MCFM, these special limits are implemented using 
a single spectator momentum, while MadDipole sums over all spectator momenta, thereby 
yielding a different subtraction term. We have checked in the relevant cases that close to 
the singularities the MCFM subtraction terms behave identical to the subtraction terms 
calculated by our code. 

We also tested the CPU time which is needed to produce the squared matrix element and 
the dipoles for a given phase space point. These checks were performed with an Intel 
Pentium 4 processor with 3.20Ghz. As an example we picked out three different processes: 

1) 99 -> 9999- \M\ 2 : 26ms, J2di P oles : 68ms 

2) uu -» ddggg: \M\ 2 : 10ms, T,dipoles : 45ms 

3) uu uuggg: \M\ 2 : 34ms, T.dipoies'- °- 15s 
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process 


subprocesses 


Drell-Yan (W) 


qq' -> W+(-> e+ Ve )g 
qg^W+(^ e+y e )q' 


Drell-Yan (Z) 


on — > Z( — > e Jr e~)Q 
qg -» e + e _ )g 


Droll- Yan (Z+jet) 


99 -» e+e-)q'q' 
qq — > Z( — > e + e~)<7Q 
qq-> Z(-» e + e~)gg 
qg -» Z(-» e+e~)qg 
gg -> e+e-)qq 


top quark pair (tt) 


qq — » t(— > bl + vi)i(—> bl~Vi)q 
qg — > t(— > bl + n)i{-^ bl~n)q 
gg -> bl + vi)t(-^ bl~Di)g 


i-channel single top 
with massive 6-quark [63] 


gg -» 

qg' — > tbq'q" 
qq' — > tbq'q" 
qg -> i6g'g 



Table 1: Set of processes for which the MadDipole code has been tested against MCFM for random 
points in phase space. All the possible initial-initial, initial-final and final-initial, dipoles for massless 
and massive final state particles have been checked with this set of subprocesses. No inconsistencies 
were found. 

The time which is needed to produce the Fortran code is strongly dependent on the process 
and ranges from a few seconds to at most a few minutes. The process 55 — > 5g is currently 
not yet feasible within MadGraph because of the size of the color factors. Once MadGraph 
has been adjusted to handle this process, it will equally become accessible for MadDipole. 

4. How to use the package 

The installation and running of the MadDipole package is very similar to the usual Stand- 
Alone version of the MadGraph code. In this section, we only provide a brief description 
for MadDipole, while more information can also be found on the MadGraph wiki page, 
http : / / cp3wks05 . f ynu . ucl . ac . be/ twiki/bin/ view/Sof tware/MadDipole. 

1. Download and extract the MadDipole package, MG_ME_DIP_V4 .4.3. tar . gz, from one 
of the MadGraph websites, e.g., http://madgraph.hep.uiuc.edu/. 

2. Run make in the MadGraphll directory. 

3. Copy the Template directory into a new directory, e.g., MyProcDir to ensure that 
you always have a clean copy of the Template directory. 

4. Go to the new MyProcDir directory and specify your process in the file . /Cards/proc_card . dat. 
This is the (n + l)-particle process you require the subtraction term for. 

5. Running . /bin/new_process generates the code for the (n + l)-particle matrix el- 
ement and for all dipole terms. After running this you will find a newly generated 
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directory . /Subprocesses/PO_yourprocess (e.g., . /Subprocesses/PO_e+e-_uuxg) 
which contains all required files. 

6. For running the check program change to this directory and run make and ./check. 

The directory . /Subprocesses/PCLyourprocess contains all necessary files needed for fur- 
ther calculations. As in the usual MadGraph code the (n + l)-particle matrix element is 
included in the file matrix. f. For the dipoles there exist several files called dipol???.f 
where ??? stands for a number starting from 001. Each file contains only one dipole. Note 
that the syntax for calling the dipoles is exactly the same as calling the (n + l)-particle 
matrix element. In particular, also the dipoles have to be called with (n + 1) momenta 
rather than with only n momenta. The remapping of the momenta going from the (n + 1)- 
particle phase space to a n-parton phase space is done within the code for the dipoles. 
The file dipolsum.f calculates the sum over all dipoles for a given (n + l)-particle phase 
space point. It contains two subroutines called DIPOLSUM and DIPOLSUMFINITE. As dis- 
cussed in Section [2,2| above, the subroutine DIPOLSUMFINITE contains the dipoles which 
only contribute potentially large logarithms but not a real singularity. 
In both subroutines the contribution of a certain dipole is only added to the sum if the 
phase space restriction specified by the a parameter is fulfilled. The value of the four a 
parameters can be changed in these two subroutines. They are all set to unity by default. 

5. Conclusions 

In this paper we have presented MadDipole, an implementation to fully automatize the 
calculation of the dipole subtraction formalism for massless and massive partons in the 
MadGraph/MadEvent framework. The implementation is done in such a way that the 
user only needs to specify the desired (n + l)-particle process and our code returns a For- 
tran code for all dipoles combined with possible Born processes which can lead to the 
(n + 1) process specified by the user. 

For the calculation of the new color factors we have used as far as possible the routines 
already provided by the original MadGraph code. We inserted the two additional operators 
for emitter and spectator and modified the evaluation of the squared color factors. 
For the contributions that are not diagonal in helicity space we again used the already 
available routines for calculating amplitudes for a given helicity combination and combined 
amplitudes with different helicity combinations to yield the off-diagonal helicity contribu- 
tions to the subtraction terms. 

We have validated the code on numerous different processes with massive and massless 
partons using two checking procedures. The ratio of A4^ +1 / ^(dipoles) was confirmed to 
approach unity as one approaches any soft or collinear limit. The package includes a file 
check_dip.f which allows the user to reproduce this check for any process under consid- 
eration. 

Secondly, we compared our code against the results of MCFM [10,11], where subtraction is 
performed using the dipole formalism, finding point-wise agreement wherever anticipated. 
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Differences with MCFM are understood to be due to different details in the implementa- 
tion. 

The MadDipole package allows the automated computation of real radiation dipole subtrac- 
tion terms required for NLO calculations. Together with the fastly developing automation 
of one-loop calculations of multi-leg processes, it could lead to a full automation of NLO 
calculations for collider processes. In view of the large number of potentially relevant 
multi-particle production processes at LHC, such automation will be crucial for precision 
phenomenology, in order to establish and interpret potential deviations from standard 
model expectations. 
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